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ABSTRACT 


Data Collection 5 processing for the Moderate Resolution Imaging Spectroradiome- 
ter (MODIS) onboard the NASA Earth Observing System EOS Terra and Aqua space- 
craft includes an algorithm for detecting multilayered clouds in daytime. The main objec- 
tive of this algorithm is to detect multilayered cloud scenes, specifically optically thin ice 
cloud overlying a lower-level water cloud, that presents difficulties for retrieving cloud 
effective radius using single layer plane-parallel cloud models. The algorithm uses the 
MODIS 0.94 pm water vapor band along with CO 2 bands to obtain two above-cloud pre- 
cipitable water retrievals, the difference of which, in conjunction with additional tests, 
provides a map of where multilayered clouds might potentially exist. The presence of a 
multilayered cloud results in a large difference in retrievals of above-cloud properties be- 
tween the CO 2 and the 0.94 pm methods. In this paper the MODIS multilayered cloud 
algorithm is described, results of using the algorithm over example scenes are shown, and 
global statistics for multilayered clouds as observed by MODIS are discussed. A theoreti- 
cal study of the algorithm behavior for simulated multilayered clouds is also given. Re- 
sults are compared to two other comparable passive imager methods. A set of standard 
cloudy atmospheric profiles developed during the course of this investigation is also pre- 
sented. The results lead to the conclusion that the MODIS multilayer cloud detection al- 
gorithm has some skill in identifying multilayered clouds with different thermodynamic 
phases. 




1. Introduction 


Plane-parallel single-layered cloud radiative transfer (RT) models are used by global 
passive imager algorithms like MODIS (MODerate resolution Imaging Spectroradiome- 
ter) (Barnes et al. 1998) for cloud thermodynamic phase, cloud-top pressure/temperature, 
and optical and microphysical properties retrievals (King et al. 2003; Platnick et al. 
2003). The use of such a RT model works reasonably well as confirmed by many field 
campaigns and theoretical calculations (King et al. 2004; Mace et al. 2005; Chiriaco et al. 
2007; Bedka et al. 2007; Otkin et al. 2008). The model can work for some retrievals if 
there are multilayered clouds in a vertical column (e.g., an ice cloud overlapping a liquid 
water cloud) and the uppermost layer is optically thick. In particular, use of the RT model 
can result in biases with cloud effective radius retrievals when liquid water clouds are 
overlaid by relatively thin cirrus clouds (Davis et al. 2009). The retrieved effective radius 
of what is thought to be single layer ice clouds decreases significantly in areas overlying 
the water clouds. When the cirrus is too optically thin to dominate the upwelling radiance 
and the cloud is identified as being liquid water phase, the retrieval tends towards abnor- 
mally large water droplets. There is not a large detrimental effect on cloud optical thick- 
ness to the extent that the combined optical thickness of all layers is retrieved with little 
dependence on the assumed phase. 

It is important to flag areas where there are problematic effective radius retrievals 
due to multilayer clouds of differing thermodynamic phases since those retrievals can ad- 
versely affect cloud statistics and should be excluded from further analysis. 

There have been other algorithms designed to identify multilayer clouds with passive 
imagers. The algorithm developed by Pavolonis and Heidinger (2004) is a pixel-level al- 
gorithm that uses ratios and differences of reflectances and brightness temperatures in 
various bands. This approach can be applied to historical and current multispectral 
imager data such as the AVHRR (Advanced Very High Resolution Radiometer) on the 
NOAA spacecraft and MODIS. Such an approach may also be continued with future 



JOURNAL OF APPLIED METEOROLOGY AND CLIMATOLOGY 


2 


measurements from the VIIRS (Visible/Infrared Imager Radiometer Suite) that will be 
flown on the NPOESS (National Polar-orbiting Operational Environmental Satellite Sys- 
tem) platforms. This multispectral approach uses ratios and differences of reflectances 
and brightness temperatures. We will show a comparison of this AVHRR- VIIRS algo- 
rithm with the MODIS algorithm in section 5. 

The algorithm developed by Baum and Nasiri (Baum et al. 2000; Nasiri and Baum 
2004) is a statistically based algorithm that is executed in shifting steps over a box area of 
user-defined size, typically 200 x 200 pixels, with a restriction that some clear sky is 
available in the area; the algorithm retrieves a probability that the cloud is multilayered. 
This algorithm was developed for the MODIS instrument, but has not been used exten- 
sively outside of case studies. The need to use a large area to work on and a requirement 
for presence of clear sky pixels within each work area reduces the effective algorithm 
resolution and usefulness as many multilayered cloud retrievals occur within synoptic 
systems that span a wide area with extensive cloud cover. We will show a comparison of 
our algorithm with the Nasiri-Baum algorithm in section 5. 

Another approach for multilayer cloud detection has been presented by Chang and Li 
(Chang and Li 2005a, b). The method of Chang and Li uses an estimation of cirrus cloud 
emissivity based on the difference of cloud top temperature retrieved by using the CO 2 
slicing result (not assuming an opaque cloud) and the 11 -pm band (i.e., assuming an 
opaque cloud). The algorithm relies on being able to identify single-layer liquid water 
clouds and clear sky pixels in an area of 250 x 250 km centered on the point of interest. 
The cloud effective emissivity is then computed, from which the infrared (IR) cloud opti- 
cal thickness is derived. If that cloud optical thickness is significantly different from the 
cloud optical thickness retrieved using a visible or shortwave infrared (SWIR) band, the 
cirrus cloud likely has a liquid water cloud underneath it. 

The MODIS operational multilayer cloud detection algorithm relies on a difference 
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in above-cloud precipitable water retrievals obtained from using the 0.94 pm band versus 
precipitable water computed from the CO 2 slicing-derived cloud top altitude. The 0.94 
pm band is relatively insensitive to optically thin cirrus and so the column moisture is 
integrated from the top of the atmosphere (TOA) to the lower level cloud, if such is pre- 
sent. The CO 2 slicing retrieval of cloud top height, and subsequent calculation of the 
above-cloud precipitable water from a forecast model profile, occurs from the TOA to the 
level of the higher cloud. From that difference, and several other tests such as the differ- 
ence between retrieved IR and SWIR cloud thermodynamic phases and reflectance ratios 
to screen for single layer clouds over bright surfaces, a determination is made as to 
whether or not the cloud is multilayered in a way that affects the applicability of the 
plane-parallel single layer cloud models used in retrievals of cloud effective radius. 

Initial CloudSat evaluations of the MODIS multilayer cloud detection algorithm 
have been done by Joiner et al.(2010) as part of a study which developed a global multi- 
layer cloud detection algorithm via cloud top pressure derived from the Aura Ozone 
Monitoring Instrument (OMI) . 

In the following discussion we present the MODIS operational multilayer cloud de- 
tection algorithm, describe how the multilayer cloud information is stored in the 
MODO 6 /MYDO 6 Level-2 HDF (Hierarchical Data Format) files, present results of exe- 
cuting the algorithm on data produced by forward simulations of multilayered clouds, and 
compare the algorithm to other methods. 

A useful cloud simulation data set was developed using a set of moist atmospheres 
from the ECMWF (European Centre for Medium-Range Weather Forecasts) 40-year 
reanalysis data set. Selected profiles were chosen at grid points that contained sufficient 
amounts of cloud to create a more realistic setting into which well-separated cloud layers 
were inserted. The profiles created in this fashion are available at http://modis- 
atmos.gsfc.nasa.gov/MOD06_L2/validation.html. 
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To summarize briefly the discussion that will follow, section 2 describes the algo- 
rithm and the data format in which the results are stored. Section 3 presents the details of 
the RT simulations and describes in detail the method used to create the simulation data 
set. Section 4 provides results from applying the MODIS operational multilayer cloud 
detection algorithm over the simulated scenes as well as selected MODIS data granules, 
and also provides an example of global statistical aggregation of multilayer cloud data. 
Section 5 discusses a direct comparison of our results with other passive remote sensor 
methods for detecting multilayer clouds. Conclusions, ongoing work and future directions 
are discussed in section 6. 

2. Algorithm description 

The operational MODIS multilayer cloud retrieval uses a number of bands in addi- 
tion to individual retrievals of physical quantities such as above-cloud precipitable water 
and cloud optical thickness, to arrive at a decision. The main component of the retrieval 
is a test for the difference of above-cloud precipitable water retrievals obtained by two 
different methods. 

The first method is based on the cloud top pressure retrieval obtained from CO 2 slic- 
ing using ratios of MODIS bands 33, 34, 35 and 36 that are centered between 13.3 and 
14.2 pm (Menzel et al. 2008). The retrieved cloud top pressure is then used to obtain 
above-cloud water vapor amount (PWC 02 ) by adding up the layer averaged water vapor 
amounts from the NCEP (National Centers for Environmental Prediction) global 6-hour 
atmospheric profile product, produced at 1° resolution. Due to the nature of CO 2 absorp- 
tion, the algorithm is sensitive to high clouds of optical thickness (t c ) greater than 0.5 
(Menzel et al. 2008) when multilayer clouds are present and will return a low value of 
above-cloud precipitable water. 

The second method uses water vapor absorption in the MODIS 0.94-pm band. 
Above-cloud precipitable water is retrieved using an iterative approach. That is possible 
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because cloud reflectance is flat in the spectral range between 0.86 and 0.94 pm and the 
difference in measured cloud reflectance is due to the water vapor amount between the 
cloud and the sensor. If the visible optical thickness of thin cirrus layer is less than 6, the 
0.94 pm band is sensitive to the low clouds when multilayer clouds are present and will 
return a higher value of above-cloud precipitable water than the CO 2 slicing method 
would. The discrepancy in retrieved amounts of above-cloud precipitable water can be 
attributed to the presence of multilayered clouds. 

The MODIS operational multilayer algorithm first assumes that a single layered 
cloud exists, with a cloud top temperature based on the 11 -pm brightness temperature. 
Cloud-top pressure is then inferred by mapping the temperature into the NCEP pressure 
profile. The mapping is done from the top downward so as to avoid the high likelihood of 
temperature inversions nearer the surface. 

This cloud top pressure together with the view geometry is used to index a MODIS 
atmospheric transmittance table for 0.86 pm and 0.94 pm, which is generated by using 
the ECMWF ERA-40 atmospheric profile database as input to MODTRAN (MODerate 
resolution atmospheric TRANsmission) version 4.2rl (Berk et al. 1998). This lookup re- 
sults in a vector of two-way atmospheric transmittance as a function of above-cloud pre- 
cipitable water for each band. These transmittance vectors are then applied to the meas- 
ured reflectances. The dominant contributor to absorption in the 0.94 pm band is water 
vapor. If there were no water vapor between observer and cloud, measured reflectances 
can be assumed to be identical. Using that assumption we look for a point where the two 
vectors intersect. The closest table index value of precipitable water at the intersection 
point is our retrieval of above-cloud precipitable water (PW0.94). We choose to neglect a 
very small amount of ozone absorption in the 0.86 pm band (<0.001 additional absorp- 
tion amount) as it has no discernible impact on location of the intersection point due to 
lookup table resolution. 
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We then use the retrieved water vapor amount to perform a crude atmospheric emis- 
sion correction on the 11 pm radiance. Measured 11pm radiance consists of three com- 
ponents: emission from ground, emission from cloud and emission from atmosphere 
above cloud. We assume that cloud emissivity is unity, therefore we do not deal with 
emission from ground. This is the exact assumption made by MODIS CO 2 slicing-based 
cloud top properties retrieval method. Now we must subtract the atmospheric emission 
from measurement and also correct the result for water vapor absorption in the 11 pm 
channel. So the final corrected radiance takes on the following form: 

, _ ~ ~ trans) 

trans 

where I meas is the measured radiance, T mean _above_cioud is the integrated layer mean 
temperature from given atmospheric profile and trans is the 1-way atmospheric transmit- 
tance at 1 1 pm. 

The entire process is repeated using the corrected 11pm radiance as a source of 
cloud top temperature. We have found that one additional iteration is enough for the re- 
trieval to converge to within 0.25K, which we consider to be a sufficient degree of accu- 
racy for our purpose. This same type of retrieval, but without iteration, is also performed 
one additional time with the assumption that the cloud in question is located at 900 hPa 
(P W 0. 94@9oo) ■ If a high, cold cloud (T c < 265 K) with little water vapor above it is moved 
vertically in the atmosphere, its retrieved temperature and pressure stay nearly constant 
because atmospheric transmittance for amounts of water vapor less than 0.5 cm shows 
very little dependence on pressure. Moving such cloud from 200 mb down to 900 mb 
changes 1 1pm transmittance by only 0.8%, 0.94pm transmittance by 1.05% and 0.86 pm 
by 0.01%. However, this is not so for a warm, low cloud with a significant amount (>1 
cm) of water vapor above it, which is fairly typical for boundary layer clouds. For such 
cloud 1 1pm and 0.86 pm transmittances change by about the same amount as for a high 
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cloud, but the 0.94 pm transmittance changes by 8% if such cloud with 1 cm of precipi- 
table water above it is moved between 600 and 900 mb. The error in retrieved precipi- 
table water amount for the lower level cloud will increase as the optical thickness of the 
overlaying ice cloud increases. The result is similar regardless of where the lower-level 
cloud lies between 800 and 1000 hPa. A low-level cloud pressure of 900 hPa is chosen as 
the default value. We mitigate the effect of ground elevation due to the fact that the 
NCEP profiles extrapolate every profile down to 1000 mb level, regardless of terrain. A 
precipitable water retrieval based on this assumption acts to mitigate the ‘cooling’ effect 
of an upper ice cloud and results in the inference of a more realistic high precipitable wa- 
ter amount above the lower level cloud. This process does not affect the results for single- 
layered ice clouds or multilayered clouds where the upper ice cloud layer is optically 
thick, and permits the tracking of more multilayered clouds. 

As both 0.94 pm and 0.86 pm channels are much more sensitive to the presence of 
lower-level clouds in multilayer situations, the retrieved precipitable water value is quite 
different from the same retrieval performed based on the inference of high clouds from 
CO 2 slicing. That difference, weighted by the total column precipitable water (TPW), is a 
key determinant of whether or not there may be multilayered clouds present. A value of 

I pw - pw I 

' _ 0.94 CQ2 ' > g Cfo 

TPW 

is used as the threshold for marking the pixel as potentially containing multilayered 
clouds based on case studies and estimates regarding the occurrence of effective radius 
biases (see following example). Forward radiative transfer simulations, discussed in sec- 
tions 3 and 4, confirm that this an appropriate choice. 

Figure 1 illustrates the effect of this retrieval on a portion of a MODIS data granule. 
These data were collected from Terra MODIS on 25 October 2008 at 0015 UTC in the 
western Pacific Ocean just east of Japan. The panels show the process of obtaining a mul- 
tilayer result using the precipitable water method. Figure la shows a false-color image of 
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MODIS bands 6, 2, 26 (1.64, 0.86, and 1.38 pm, respectively). Thin cirrus is advecting 
over a field of cumuliform clouds. In this false-color composite, liquid water clouds ap- 
pear gold, ice clouds appear blue and white, and the ocean surface appears black. A num- 
ber of areas where thin cirrus overlaps the liquid water clouds are visible in the image and 
take on a greenish hue. Figure lb is an image of above-cloud precipitable water from the 
MODIS cloud top properties algorithm that uses CO 2 slicing (PWco 2 )- The figure shows 
a strong preference for high cloud properties in the overlap region, and thus a low water 
vapor amount is derived above clouds. There is barely a trace of the low-level clouds in 
the image. Very low values of precipitable water are seen for the high clouds as expected. 
Figure lc is an image of the standard 0.86-0.94 p m retrieval of precipitable water 
(PW0.94), which is more sensitive to low clouds and so gives higher precipitable water 
values that are more typical for those clouds. Figure Id is the difference image between 
the precipitable water from the CO 2 slicing and the 0.94 pm algorithm. Outlines of low- 
level clouds are becoming clearly visible in the difference image. Small differences in 
precipitable water correspond to either thicker cirrus, which is not sensitive to multilayer 
clouds, or breaks in the low-level cumulus clouds. But more cloud could be flagged as 
the cirrus becomes thicker to the west and is affecting the vertical placement of the cumu- 
lus. Figure le shows the precipitable water retrieval in which the low-level clouds are 
assumed to be at the 900 hPa level. It is not that different from the main 0.86-0.94 pm 
result with the exception that it captures some of the cloud features covered by somewhat 
thicker cirrus to the west. Even though the clouds are thicker, they still contain some con- 
tribution from the underlying low-level cloud. Figure If shows the difference image re- 
sulting from the 900 hPa retrieval versus that from the CO 2 slicing. 

The final two images are the retrieved cloud optical thickness and effective radius 
for the scene. The warm colors indicate liquid water clouds with cold colors for ice cloud 
retrievals. The optical thickness image indicates that the cirrus is quite thin and fairly uni- 
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form over the overlap area. There is no significant impact of multilayered clouds on opti- 
cal thickness as the overlying cirrus is thin and its contribution to the combined visible 
optical thickness is very small. In contrast, the impact on the cloud effective radius re- 
trieval is much greater. The outlines of low-level clouds are clearly seen in the effective 
radius image as areas of small ice effective radii. The breaks of open water in the cumu- 
lus cloud fields return effective radius values of around 25 pm, so it is unlikely that the 
actual cloud microphysics is changing in the overlap area. 

Figure 2 shows the net statistical effect of multilayer clouds on cloud optical thick- 
ness and cloud effective radius. While there is not a large effect on cloud optical thick- 
ness, there is a significant shift in effective radius distribution towards smaller radii when 
multilayered clouds are not removed from the scene. In this particular case 19.2% of ice 
cloud in the scene was multilayer and -54,000 pixels were removed from the distribution. 

The MODIS CO 2 slicing algorithm is applied with the most confidence for clouds at 
pressures lower than about 700 hPa (Menzel et al. 2008). In a typical MODIS scene, 
however, the CO 2 slicing algorithm is rarely applied for clouds at pressures larger than 
600 hPa. If the CO 2 slicing algorithm is unable to converge on a solution, the 1 1-pm band 
is used under the assumption that there is a low-level opaque cloud present. The choice 
was made to ignore CO 2 slicing results at pressures larger than 550 hPa to minimize the 
potential for false positive retrievals. In light of improvement in vertical resolution to 101 
levels used in MODIS CO 2 slicing algorithm beginning with Collection 6, this 550 hPa 
restriction may be eased in the future, although uncertainties due to resolution of the 
NCEP profiles will remain. 

Due to uncertainties in inferring cloud emissivity from passive sensors, it is possible 
to obtain a false positive multilayer retrieval for the case when an optically thin cirrus 
cloud is present with x c < 4. If the cloud is very optically thin, upwelling radiance from 
surface will cause that cloud to be placed at pressure much higher than truth. That means 
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the 0.94 pm cloud top properties method will retrieve much higher precipitable water 
amount than CO2 slicing would because of surface contamination and not because of 
multilayer situation. We assume that if the total column optical thickness is < 4, the like- 
lihood is that there is not a lower cloud underneath it. The liquid water cloud layer under- 
neath would most likely push the total optical thickness above 4. If a liquid water cloud is 
so thin that threshold of 4 is not reached, then we would have difficulty with retrieving 
effective radius due to shape of forward library space (Platnick, et.al. 2003), any multi- 
layer situation aside. False negatives do arise from use of this threshold, but with overall 
effective radius retrieval uncertainty being well above 20% for thin clouds, the weight of 
such retrievals should be greatly reduced in any statistical studies anyhow. 

We also must consider cases of single-layer clouds over bright surfaces. It is possible 
for the algorithm to mistake a thin cirrus cloud over a bright surface for a cloud that is 
multilayer. The 0.65 pm and the 1.24 pm reflectances are used to check for vegetation 
and snow/ice, respectively. Cloud reflectance is reasonably flat in that spectral region, 
while surface albedo changes significantly. So for a true multilayer cloud situation, the 
reflectance ratio would be close to 1.0, but not so for a single layer of thin cirrus over a 
bright surface. It is useful to use ratios of 0.86 pm reflectance to 0.65 pm and 1.24 pm 
reflectance to check for bright surfaces, with thresholds set as follows: 

/?(0.86 pm) < | 25 
R( 0.65 pm) 

R( 0.86 pm) < | j 
/?(1.24 pm) 

These thresholds were empirically derived on the basis of case studies; however our for- 
ward simulations indicate that a parameterization based on ecosystem type may be more 
appropriate in the future. We will investigate such parameterization in MODIS data for 
collection 6. 

In addition to the precipitable water difference, another test is based on retrievals of 
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cloud thermodynamic phase from two different methods. The first method is the MODIS 
SWIR thermodynamic phase (SP) algorithm (Platnick et al. 2003) that uses a number of 
cloud mask tests and reflectance ratios in visible, NIR and SWIR bands to arrive at cloud 
thermodynamic phase. The second method is the IR bi-spectral cloud phase (IP) algo- 
rithm based on brightness temperature differences between 8.5 and 1 1 pm bands, which 
is a modification of the Baum IR tri-spectral algorithm (Baum et al. 2000). When these 
two methods infer different thermodynamic phases, that can be an indication of a multi- 
layered cloud situation. This particular test tends to be sensitive to cirrus over liquid wa- 
ter clouds in which thin cirrus is too thin to result in an ice phase retrieval, but still biases 
the liquid water cloud retrievals as the cloud effective radius retrieval is larger than ex- 
pected. 

The main uncertainty associated with using the thermodynamic phase test tends to 
arise in polar regions. At latitudes above 60°, the IR method results in quite a few unde- 
termined phase answers due to inherent difficulties of an IR method over very cold sur- 
faces, so we assign a lower degree of confidence to multilayered clouds that are flagged 
only by the cloud phase test and no other test. 

The 0.94 pm precipitable water retrieval performed at both pressure at cloud top and 
at 900 mb, together with a test on retrieved cloud thermodynamic phase combine to cre- 
ate a final integer answer that tells the user whether the multilayer detection algorithm 
arrived at a positive result and what method(s) were positive as shown in Table 1. We 
store the final value in the MOD06/MYD06 Level-2 HDF file in two places. The values 
from table 1 are stored in a Scientific Data Set (SDS) named Cloud Multi Layer Flag. 
The multilayer cloud information is also stored in the 5 th byte of the Qual- 
ity Assurance lkm SDS as information about the thermodynamic phase of the cloud and 
its multilayer status. The full description of the Quality Assurance lkm SDS is given in 
Hubanks (2006) and a brief listing of relevant values is given in table 2. 
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The discussion in this section is summarized in Figure 3. The algorithm flow chart 
shows the overall logical flow of the algorithm. 

3. Radiative transfer models 

We have conducted an extensive set of forward RT modeling studies of multilayer 
clouds under varying atmospheric conditions, layer separations, surface types and layer 
thicknesses to thoroughly test the sensitivities and skill of the MODIS multilayer cloud 
detection algorithm. 

Zonal and temporal average profiles are calculated from the ECMWF sampled 60- 
level global atmospheric profile database aggregated from ERA-40 data over 48 days for 
two years using 1 st and 1 5th of each month between January 1 992 and December of 1 993 
(Chevallier, 2001). The database profiles were separated to represent a typical midlati- 
tude summer (MLS), midlatitude winter (MLW), tropical atmosphere (TRP), and polar 
oceanic (POL) profile. Profiles over polar landmasses, dominated by profiles from the 
Antarctic continent, were not included as they would contain strong inversions and would 
be likely used disproportionately for pressures lower than 700 hPa. The polar oceanic 
profile consists of daytime profiles only. Nighttime profiles are not used since for our 
purposes, cloud optical and microphysical property retrievals are performed in daytime 
only. We define the tropical region as 30°S < latitude < 30°N, midlatitudes as 30° < 
|latitude| < 60° and the polar regions as above 60° latitude. For midlatitudes, winter pro- 
files occur between 1 November and 30 April; summer profiles are the remainder of the 
year. Within each latitude belt, profiles are chosen from regions that had cloud fraction 
(CF) > 0.85 to match the conditions of interest. Profiles were separated further by land 
and ocean using the ECMWF land fraction flag with threshold set at 0.5. 

Given these averaged profiles, chosen levels were saturated with cloud having an 
appropriate cloud thermodynamic phase by setting relative humidity at levels that were 
chosen to contain cloud to 100%. The profiles were interpolated from the native 60-level 
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resolution to 36 levels spaced at 1 km vertically between 0 and 25 km with sparser reso- 
lution in the upper atmosphere. 

Figure 4 shows a combined plot of the temperature and moisture profiles used in the 
simulations. These particular plots show the liquid water cloud layer at 2 km. Simulations 
were run for a variety of solar and view zenith angles with the solar zenith angles appro- 
priate for the time of year in question. We sampled the solar zenith angle from the 
MODIS Level-3 global monthly product (Hubanks et al. 2008). The cosine of the view 
zenith angle corresponded to p = 1.0, 0.8 and 0.6. For detailed examination, simulations 
were run for ice cloud effective radii of 10, 30 and 50 pm and water radii of 6, 10 and 20 
pm. An ice cloud layer of 2 km physical thickness was fixed at the base of the tropopause 
as indicated by temperature in each of the different profiles shown in Figure 4: 8 km 
(MLW and POL), 12 km (MLS) and 14 km (TRP) . The ice cloud optical thickness var- 
ied between 0 and 20, with increments as follows: 0.0, 0.1, 0.25, 0.5, 1.0, 2.0, 3.0, 4.0, 
5.0, 6.0, 7.0, 8.0, 10.0, 20.0. These increments were chosen specifically to examine the 
thin cirrus region and also to appropriately capture the point where the ice cloud becomes 
too thick to be affected by the underlying liquid water cloud. Water cloud layers were 
assumed to be 1 km thick and were placed at two different altitudes: 2 km and 4 km. For 
liquid water clouds, optical thicknesses ranged as follows: 0.0, 2.0, 5.0, 10.0, 20.0. 

Radiances were simulated for 16 MODIS bands, which was necessary to perform the 
relevant cloud mask tests (Ackerman et al. 2006; Frey et al. 2008), in particular the 3.7- 
1 1 pm brightness temperature test, the CO 2 slicing cloud top properties retrieval, and the 
full MODIS cloud optical and microphysical property retrievals. The set included 
MODIS bands 1, 2, 5, 6, 7, 19, 20, 22, 26, 29, 31, 32, 33, 34, 35 and 36 (Ackerman et al. 
2006). 

Each simulation was repeated over a wide variety of surfaces. The oceanic profiles 
only had one option (dark ocean with surface albedo of 0.05) with the exception of polar 
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ocean that also included a sea ice surface. The land surface profiles presented options of 
vegetated, desert, or snow cover. Midlatitude land included mixed forest and desert with 
or without snow, appropriately, while tropical land included desert and evergreen broad- 
leaf forest. All classifications were based on definitions of the IGBP ecosystem map and 
the surface albedo values taken from MOD43 -based 1-km resolution surface albedo 
product (Moody et al. 2005, 2007, 2008). Figure 5 shows a plot of the white-sky (diffuse) 
surface albedo as a function of wavelength for the various surfaces considered in this in- 
vestigation. MODIS bands that contain no solar component were given a zero surface al- 
bedo. 

The RT simulations were performed using DISORT (Discrete Ordinate Radiative 
Transfer) code (Stamnes et al. 1988) using liquid water cloud phase function results from 
Mie calculations based on the water droplet size distributions using a gamma distribution 
with an effective variance of 0.1 (Platnick et al. 2003) and bulk ice cloud phase functions 
developed by Baum et al. (2005a, b). The same phase functions for both ice and liquid 
water are used in the LUTs employed in the MODIS cloud optical and microphysical 
properties algorithm for collection 5. The correlated-^: method (Kratz, 1995) was used to 
account for water vapor and other gaseous absorbers. The DISORT code, in conjunction 
with the correlated- A: method, then produced the simulated MODIS band radiances. We 
used 32 streams in our radiative transfer calculations, which, together with truncation of 
strong forward peaks and use of delta-fit method by Hu et al. (2000), can be considered 
sufficient computational accuracy as described by Ding et al. (2009). 

With the parameter ranges described above, the forward RT calculations resulted in 
26 files corresponding to combinations of atmospheric profiles and surface types. Each 
file contained 7560 individual data points for each geometry, optical thickness and 
effective radius tested. Results are provided in the following section for application of the 
MODIS multilayer cloud detection algorithm to a cross-section of this database of 
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simulated MODIS radiances. In section 5 similar results are provided for the Pavolonis- 
Heidinger and Nasiri-Baum algorithms to this same dataset with comparison to the results 
from the MODIS operational algorithm. 

4. Results 

In this section we show results of applying the MODIS multilayer cloud retrieval 
simulated MODIS data. We show the results from a cross-section of our forward RT 
simulations. Figure 6 shows a set of combined results from the DISORT forward simula- 
tions. To facilitate the interpretation of results, we group individual runs having all but 
one identical parameters to illustrate the effect of the differing parameter on the multi- 
layer cloud retrieval result. Figure 6a combines the results of simulations conducted with 
a nadir view, solar zenith at 32°, dark ocean surface and liquid water cloud located at an 
altitude of 2 km. The atmospheric profile is varied in terms of the overall column mois- 
ture content. The plot in Figure 6a effectively shows multilayer cloud detection as a func- 
tion of the total column water vapor. The ‘bits’ in the effective binary numbers that result 
from this data combination indicate whether or not a multilayer cloud was detected. The 
bit significance was arranged as a function of the column moisture with the least signifi- 
cant bit for the most moisture. For example, a value of 01 1, which is light green in the 
plot, means that a multilayer cloud was detected under the conditions specified above us- 
ing TRP and MLS profiles, but no multilayer cloud was detected for the MLW profile. 
The algorithm is more likely to detect a multilayer situation when the ice cloud is opti- 
cally thin if the atmospheric moisture content is higher. 

Figure 6b shows the same basic situation as Figure 6a with the exception that the al- 
titude of the lower-layer liquid water cloud was placed at 4 km and thus decreases the 
cloud layer separation. When the cloud layer separation is smaller, the amount of atmos- 
pheric water vapor between the two cloud layers is also lower and so the absorption in the 
0.94 pm channel is decreased over the previous case. The sensitivity of the algorithm de- 



JOURNAL OF APPLIED METEOROLOGY AND CLIMATOLOGY 


16 


creases as the ice optical thickness increases compared to the case where the liquid water 
cloud is at 2 km altitude. Some false positives occur in which multilayer cloud is detected 
for thicker liquid water clouds where there is no ice cloud above. These false positives 
come from growing uncertainties in retrieving IR cloud phase and CO 2 cloud top proper- 
ties as the cloud gets colder. The detection results can be inspected further by looking at 
individual tests, some of which have lower confidence than others as mentioned in sec- 
tion 2. The detection status is reported as a binary answer and may result in a false posi- 
tive result. 

Figure 6c illustrates the multilayer detection result as a function of underlying sur- 
face type, assuming a single MLS profile and a liquid water cloud placed at 2 km altitude. 
The surface types are arranged such that the least significant bit corresponds to the lowest 
overall surface albedo with no snow on the ground. The plot shows that multilayered 
clouds are not detected for a desert ecosystem with thin liquid water clouds below, since 
the liquid cloud emissivity is likely somewhat less than 1 .0, thereby indicating that we 
may need a separate detection threshold for deserts since the surface albedo of deserts is 
significantly different in spectral shape from vegetation and snow/ice surfaces. The desert 
spectral albedo tends to be somewhat flatter than vegetation, as Figure 5 shows, and so 
may require a somewhat different approach. The effect of this on our global statistics is 
not very significant as the actual cloud fraction over deserts is rather low (cf. Figure 9). 
Figure 6d shows multilayered cloud detection as a function of cosine of the viewing ze- 
nith angle (p) for a MLS profile with a dark ocean surface and a lower-layer liquid water 
cloud placed at 2 km. The points are ordered in p -space such that a more oblique angle, 
i.e., lower p, is the least significant bit in the binary number displayed. The relative azi- 
muth angle for this comparison was set to 0°. The figure indicates that the algorithm is 
more likely to detect a thinner ice cloud over a liquid water cloud at more oblique angles. 
On the other hand, it is possible to flag cases with higher ice cloud optical thicknesses at 
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more nadir view angles. 

The Pavolonis-Heidinger method, originally developed for the AVHRR and adapted 
for the upcoming VIIRS instrument, uses a series of reflectance and brightness tempera- 
ture difference thresholds described in detail in (Pavolonis et al. 2004). For the algorithm 
comparison purposes we have been provided with their most recent development of the 
method, with improvements and modifications made since the publication of their paper. 
Similar to the MODIS operational algorithm, it is a single-pixel method that works on 
samples individually without using any spatial aggregation. Because of this similarity we 
were able to execute the Pavolonis-Heidinger algorithm on the results of our DISORT 
simulations of multilayer clouds. Figure 7 shows the comparison of these results. The 
figures on the left are the MODIS results from Figure 6a, 6b and 6d, and on the right are 
corresponding Pavolonis-Heidinger results. We compared the algorithms for three out of 
four database cross-sections shown in Figure 6. It was not possible to perform the exact 
comparison for the surface-type section, since the Pavolonis-Heidinger algorithm uses a 
0.41 pm band over desert regions that we did not include in the original DISORT band 
set. The Pavolonis-Heidinger algorithm uses a lookup table (LUT) derived from simula- 
tions of multilayered clouds over a various surfaces. The LUT includes the difference in 
brightness temperatures (BTD) between the 11 and 12 pm bands. A threshold function is 
defined since the multilayered clouds (i.e., ice over water cloud) display a BTD as a func- 
tion of visible reflectance that is quite different from single-layered liquid water and ice 
clouds. In addition to that threshold, a number of constraints are placed on reflectances at 
0.65 and 1.38 pm to help with the identification of single layer clouds over a variety of 
surfaces. The 1 .65 pm band is used by the algorithm to aid in identifying the thermody- 
namic phase of clouds since ice clouds have greater absorption than liquid water clouds at 
1.65 pm. 

There are similarities in the results as well as some differences, but overall the com- 
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parison is favorable. The MODIS algorithm has a somewhat wider section where multi- 
layer clouds are detected for the entire range of the varied conditions, be it atmospheric 
moisture content or view angle. However the detection rate generally drops off as the ice 
cloud thickens with only the thickest simulated liquid water cloud showing at ice cloud 
optical thickness of 10. Both algorithms show that once the ice cloud optical thickness 
reaches 20, no detection of multilayer clouds is possible. Both algorithms also show that 
detection is a function of layer separation with detection rate being lower when the liquid 
water cloud is placed at 4 km as opposed to 2 km cloud top altitude. The Pavolonis- 
Heidinger algorithm shows more detection when both cloud layers thicken, but not as 
much when the cloud layers are thin. 

We could not run the Nasiri-Baum algorithm on DISORT simulations because it is a 
statistical aggregate algorithm that depends on natural variability of the data within a 
given area. More specifically, the algorithm uses a brightness temperature difference be- 
tween 8.5 and 11 pm to confidently determine the clouds that are liquid water phase. 
Liquid water clouds tend to have a large, negative brightness temperature difference. The 
data are analyzed on 200 x 200 pixel tiles that must meet a number of conditions to at- 
tempt the retrieval. The tile must contain clouds and must contain at least 10 pixels each 
of ice cloud, liquid water cloud, and clear sky. Additionally, a scatter plot of 11 -pm 
brightness temperature versus the 2.1 pm reflectance is created, wherein pixels that be- 
long to single layer ice clouds and single layer liquid water clouds create two distinct 
lines. The clear sky points lay on the intersection of those lines. The angle between the 
lines must be greater than 20°. The grid of tiles undergoes successive systematic shifts to 
increase the number of times a particular pixel is processed; the more times a pixel is 
flagged as multilayered, the higher the confidence of the final answer. As DISORT re- 
sults are single points, there is no way to create an appropriate analysis box that would 
satisfy the data requirements of the Nasiri-Baum algorithm. 
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Our overall conclusion from examining all these results is that the MODIS multi- 
layer cloud detection algorithm is robust and performs as intended under a wide variety 
of conditions. 

5. Analysis and comparison with other methods 

In this section we show an example case study from MODIS and comparisons of our 
method against two other multilayer cloud detection algorithms, which we mentioned in 
section 1. 

Figure 8 shows an example of multilayer cloud detection for a Terra MODIS granule 
acquired on 25 October 2008 off the coast of Japan at 0015 UTC. This is a full granule, a 
portion of which was shown in Figure 1 . Figure 8a shows an atmospherically corrected 
true color image formed as a composite of MODIS bands 1, 4, and 3 (0.65, 0.55, and 0.47 
pm, respectively). While the false color image indicates where the clouds are, it provides 
very little information about the various cloud layers or the thermodynamic phase of the 
clouds. Figure 8b shows a false color image formed as a composite of bands 6, 2, and 26 
(1.64, 0.86, and 1.38 pm, respectively), which more readily separates clouds of different 
thermodynamic phase by color. There is a significant amount of multilayer cloud in this 
scene, indicated by areas where the yellow liquid water clouds show through the more 
blue and white ice clouds. Figure 8c shows the results from applying the multilayer cloud 
detection algorithm. Different values on the color scale correspond to tests flagging the 
cloud as clear sky (0), single layer cloud (1), and multilayer (2-8) cloud, as described in 
Table 1. These results are not an absolute measure of multilayer cloud amount, but rather 
provide a map of areas where the presence of multilayer clouds adversely affects cloud 
effective radius retrievals. 

In Collection 5, MODIS multilayer cloud retrievals are aggregated to the global 
level-3 daily, eight-day, and monthly products as an average of data down-sampled to 5 
km and aggregated into a 1 ° grid. The multilayer cloud fraction is stored, combined and 
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separated by thermodynamic phase, and also includes mean values of cloud optical and 
microphysical properties retrievals, both with and without multilayer clouds. (Hubanks et 
al. 2008). Figure 9 shows an example of such an aggregation for the month of October 
2008 derived from Terra MODIS data. Figure 9a shows the fraction of all cloudy pixels 
that have the multilayer flag set, and Figure 9b shows the mean monthly cloud fraction. 
The small black area on the very top of the images corresponds to polar darkness or low 
sun where no retrievals are attempted (cosine of the solar zenith angle po < 0.15). 

A monthly global map like this is useful for providing the spatial distribution of mul- 
tilayered clouds. Based on observational evidence, one might expect a higher frequency 
of multilayered clouds to occur in the vicinity of low-pressure systems and their frontal 
boundaries. Higher frequencies of multilayer clouds tend to occur in the Southern Ocean 
and in the storm tracks (higher latitude zones) of both hemispheres. Tropical anvil cirrus 
is also a likely candidate to create multilayer cloud situations. A good portion of the In- 
tertropical convergence zone (ITCZ) is flagged as multilayer in the eastern Pacific Ocean. 
Strong convective zones over rainforest areas also tend to generate anvil cirrus, resulting 
in high frequencies of multilayered clouds in the Congo basin, Borneo, and New Guinea. 
One can also note the effect of advection of anvil cirrus over the marine stratocumulus 
zones in the Southern Hemisphere off the coasts of Peru and Ecuador, and in the Gulf of 
Guinea. 

The three multilayer cloud detection algorithms previously discussed are now ap- 
plied to the MODIS granule shown in Fig. 8, with the results shown in Figure 10. Figure 
10a shows the true color composite constructed from bands at 0.65, 0.55, and 0.47 pm, 
Figure 10b the false color composite constructed from bands at 1.64, 0.86, and 1.38 pm, 
and Figure 10c the false color composite constructed from bands at 0.55, 1.64, and 2.13 
pm. Figures lOd-f show the results of applying the multilayer cloud detection using the d) 
MODIS operational algorithm, e) Pavolonis-Heidinger algorithm, and f) Nasiri-Baum 
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algorithm. 

As there is a wide range of options, described in the code documentation, that the 
Nasiri-Baum algorithm can be executed under, for the purposes of this comparison we 
took the suggested default values. The Nasiri-Baum algorithm can only be executed un- 
der conditions that some clear sky, liquid water cloud and ice cloud exists within the box 
being currently analyzed, so the algorithm does not attempt retrievals over a portion of 
this granule. The Nasiri-Baum algorithm also outputs its result as a probability of multi- 
layer cloud being present. For clarity we display only a non-zero overlap probability as a 
positive answer. We performed a similar procedure with the results from the MODIS op- 
erational multilayer cloud algorithm, combining the multilayer values 2-8 into a single 
positive identification value. The Pavolonis-Heidinger algorithm returns its result as a 
single value so no additional data conversion was necessary to visualize the results. 

Overall many of the same areas flagged as multilayer, even though the results may 
not look exactly the same, as the different multilayer algorithms were developed with dif- 
ferent applications in mind. The Nasiri-Baum algorithm gives the fewest multilayer oc- 
currences, but that can be attributed to a limited area over which the algorithm attempts 
retrievals. The main disagreement between Pavolonis-Heidinger and our algorithm arises 
in the flagging of thicker high clouds as being part of multilayer scenes (e.g., left side 
portion of the granule). One might argue that most, if not all, ice phase clouds in that part 
of the granule are multilayered clouds because of the apparent wide presence of low 
clouds in that region as well as there being some indication in the 1.38 pm false color 
composite. The result given by the Pavolonis-Heidinger algorithm is consistent with de- 
tection achieved for simulated DISORT data, where clouds with combined extinction op- 
tical thickness as large as 30, with the upper layer thickness of 10, can be flagged as mul- 
tilayer, as shown in Figure 7. 

The decision whether to flag a cloud as multilayer depends on the issue being ad- 
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dressed. In our case, we are looking for multilayer clouds that challenge the applicability 
of our single-layer plane-parallel cloud models used in cloud optical and microphysical 
property retrievals. Our goal is to create a map of areas where the model application is 
problematic. From our RT simulations we have found that the effect of ice cloud over- 
lapping a liquid water cloud on cloud effective radius retrieval diminishes quite rapidly 
with increasing ice cloud optical thickness and is barely detectable when ice cloud optical 
thickness becomes greater than about 6. While it may be the case that the thicker upper 
level clouds in this granule are also multilayer, having those clouds flagged as such does 
not impact our primary objective regarding microphysical biases. 

Joiner et al. have provided initial comparisons of the MODIS multilayer cloud detec- 
tion algorithm to CloudSat/CALIPSO as part of a study which developed a multilayer 
detection algorithm via OMI-derived cloud-top pressure (UV rotational-Raman scatter- 
ing) compared with MODIS thermal emission retrievals. The OMI approach is philoso- 
phically similar to the MODIS approach reported here in that both take advantage of solar 
reflectance path lengths being affected by gas species between cloud layers (O 2 for OMI, 
water vapor in this study) in a fundamentally different way than CO 2 slicing spectral 
bands; an O 2 method of course have the advantage of being a well-mixed gas (see Fig. 6a, 
b). With multilayer detection from CloudSat defined as layer separations greater than 
200 hPa, global analyses from a single day in Joiner et al. found that the MODIS algo- 
rithm correctly identified single and multilayer cloud layers 83.4% of the time, with false 
positives and negatives, 9.8% and 6.8% of the time, respectively. However, the MODIS 
algorithm is sensitive to upper layer cirrus optical thicknesses on the order of 0.2 (Figs 6, 
7). False positives from missed thin cirrus are not considered in the Joiner et al. study. As 
already noted, the MODIS algorithm was intended to flag cases where upper layer clouds 
have optical thicknesses too small to significantly screen lower level clouds in the short 
wave infrared bands used for particle size retrievals; such cases were not considered in 
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assessing false negatives. For both types of false detection, information from CALIPSO 
cirrus optical thickness retrievals are needed, as well as monthly and seasonal compari- 
sons (vs. single day). 

6. Conclusions and future directions 

In this paper we present the MODIS operational multilayer cloud detection algorithm 
used in the MODIS collection 5 cloud optical and microphysical properties product 
(MOD06 for Terra, MYD06 for Aqua). The multilayer cloud detection method was de- 
veloped to address a need to indicate areas of cloud where an assumption of single-layer 
plane-parallel cloud models was challenged due to the presence of two distinct cloud lay- 
ers with differing thermodynamic phases, with the upper cloud layer being optically thin. 
Such situations manifest themselves as areas of abnormal cloud effective radius retriev- 
als. Our method uses the difference between retrieved above-cloud precipitable water 
amounts from the 0.94 pm band and from the CO 2 slicing cloud top height, together with 
a number of other tests. The physical basis of the multilayered cloud detection algorithm 
is provided, with examples of results from forward simulations as well as case studies 
involving MODIS data and global aggregations of results. Results from this approach are 
compared to two other methods of multilayer cloud detection. We also present a set of 
standard cloudy atmospheres that we developed to perform our studies. Wherever possi- 
ble we perform all comparisons using a single source dataset, so the differences in re- 
trieved results are solely due to differences in methodology. 

Our results and analysis indicate that the multilayer cloud detection algorithm pre- 
sents a reliable means of identifying situations that would create difficulties for retrievals 
of cloud effective radius. The forward simulations indicate that there are very few false- 
positive results and that they arise under conditions that would result in high retrieval un- 
certainty due to one of the cloud layers being extremely thin. Forward radiative transfer 
simulations, performed under a wide variety of surface and atmospheric conditions, are 
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used in our analyses to provide further insight as to the robustness of the algorithm. 

We are currently investigating a number of improvements for the MODIS opera- 
tional multilayer cloud detection algorithm that may be implemented for MODIS collec- 
tion 6. Those improvements involve bringing in additional retrievals of physical quanti- 
ties performed using different methods, which in our experience has shown to contain 
multilayer cloud information. We also intend to continue our ongoing comparisons by 
performing more extensive MODIS to CALIPSO comparisons. However for that work 
we require more data products than what is currently available from CALIPSO. We have 
begun such studies (R. Holz, private communication) and are awaiting the next release of 
CALIPSO products (Version 3). 
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TABLE LEGENDS 

Table 1. Listing of discrete values in SDS Cloud Multi Layer Flag and definitions. 
Table 2. Listing of discrete values in the 5 th byte of SDS Quality Assurance lkm and 
definitions. 
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Table 1. Listing of discrete values in SDS Cloud Multi Layer Flag and definitions. 


Result 

Value 

Description 

0 

Clear sky 

1 

Single layer cloud or cloud too thin (x c < 4) 

2 

Multi-layer. Cloud phase test positive 

3 

Multi-layer. Precipitable water with retrieved pressure test positive 

4 

Multi-layer. Precipitable water with pressure fixed at 900 hPa test positive 

5 

Multi-layer. Both 3 and 4 

6 

Multi-layer. Both 2 and 3 

7 

Multi-layer. Both 2 and 4 

8 

Multi-layer. All three tests positive 
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Table 2. Listing of discrete values in the 5 th byte of SDS Quality Assurance lkm and 


definitions. 


Result 

Value 

Description 

0 

Cloud mask undetermined 

1 

Not processed (typically clear) 

2 

Single layer liquid water cloud 

3 

Multilayer liquid water cloud 

4 

Single layer ice cloud 

5 

Multilayer ice cloud 

6 

Single layer undetermined phase cloud 

7 

Multilayer undetermined phase cloud 
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Fig. 1. 


Fig. 2. 


Fig. 3. 
Fig. 4. 


FIGURE LEGENDS 

Cloud and water vapor properties over the western Pacific Ocean off Japan as 
acquired by Terra MODIS on 25 October 2008 at 0015 UTC. Panel (a) is a false 
color composite of one MODIS granule, showing liquid water clouds in gold, 
ice clouds in blue and white, and overlapped clouds in green. Panels (b) and (c) 
show above cloud precipitable water derived from CO 2 slicing and 0.94 pm so- 
lar absorption, respectively. Panel (d) shows the difference in the derived PW 
from these two techniques. Panels (e) shows above cloud PW derived by as- 
suming the reflecting lower cloud is located at 900 hPa and (f) the difference in 
above-cloud PW from this result and CO 2 slicing. Panels (g) and (h) show re- 
trieved values of cloud optical thickness and effective radius for liquid water 
and ice clouds. Panel (i) shows results of the precipitable water multilayer tests. 
Areas of color indicate where the tests returned a positive answer. I.e. precipi- 
table water difference was higher than 8% of the total column water vapor. 
Histograms of optical thickness and effective radius for ice clouds within the 
scene presented in Figure 1 . The effective radius histogram shows a significant 
‘shoulder’ of smaller effective radii when multilayer clouds are not removed 
from the scene. Those small effective radius retrievals come from the mixing of 
the strong liquid water cloud signal with a relatively weak cirrus signal. 
Flowchart for determining the presence of multilayer clouds using MODIS (col- 
lection 5). 

Vertical profiles of (a) temperature and (b) moisture used as a base in the for- 
ward models. The moisture profiles show saturation at levels that contain 
clouds. In this example the liquid water cloud is located at 2 km altitude and ice 
clouds are placed according to the tropopause location. 

White-sky albedo as a function of wavelength for selected IGBP ecosystem 
classifications used in the forward calculations. 


Fig. 5. 
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Fig. 6. MODIS multilayer cloud detection over various surfaces, water vapor content, 
and view zenith angle for a cross-section of DISORT simulations. . Water vapor 
profiles are mid-latitude summer (MLS), mid-latitude winter (MLW) and mean 
tropical (TRP). Panel (a) indicates the detection of multilayer clouds over ocean 
when the lower layer liquid water cloud is placed at 2 km altitude and panel (b) 
for a lower layer liquid water cloud placed at 4 km altitude. The color bar iden- 
tifies the selection of multilayer clouds for various atmospheric profiles contain- 
ing different water vapor profiles. Panel (c) shows the MODIS multilayer cloud 
detection over various surfaces including desert, forest, and snow. Panel (d) 
shows multilayer cloud detection for various view zenith directions. 

Fig. 7. Results of MODIS (left) and Pavolonis-Heidinger (right) multilayer cloud de- 
tection for a cross-section of DISORT simulations. Panel (a) indicates the detec- 
tion of multilayer clouds over ocean when the lower layer liquid water cloud is 
placed at 2 km altitude and panel (b) for a lower layer liquid water cloud placed 
at 4 km altitude. The color bar identifies the selection of multilayer clouds for 
various atmospheric profiles containing different water vapor profiles. Panel (c) 
shows multilayer cloud detection for various view zenith directions. 

Fig. 8. Multilayer cloud over the western Pacific Ocean off Japan on 25 October 2008. 
Panel (a) is a true color composite of a Terra MODIS granule at 0015 UTC. 
Gray box indicates the area of interest presented in Figure 1. Panel (b) is a false 
color composite that more clearly shows the liquid water clouds (gold), ice 
clouds (blue) and multilayer clouds (green). Panel (c) is the multilayer cloud 
SDS that shows clear sky (0), single layer clouds (1), and multilayer clouds de- 
termined by various algorithm choices (2-8) (see Table 1). 

Fig. 9. Terra MODIS monthly level-3 global product for October 2008. Panel (a) 
shows the fraction of cloudy pixels with multilayer cloud flag identified (all 
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phases), and panel (b) shows the combined cloud fraction. 

Fig. 10. Multilayer cloud analysis and cloud optical properties over the western Pacific 
Ocean off Japan as acquired by Terra MODIS on 25 October 2008 at 0015 
UTC. Panel (a) is a true color composite of one MODIS granule and panel (b) is 
a false color composite, showing liquid water clouds in gold, ice clouds in blue 
and white, and overlapped clouds in green. Panel (c) is another false color com- 
posite showing ice clouds in pink and liquid clouds in blue. Panels (d), (e), and 
(f) show multilayer cloud identification using the MODIS collection 5, Pavolo- 
nis-Heidinger, and Nasiri-Baum algorithms, respectively. Panels (g) and (h) 
show retrieved values of cloud optical thickness and effective radius for liquid 
water and ice clouds. 
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Fig. 1 . Cloud and water vapor properties over the western Pacific Ocean off Japan as 
acquired by Terra MODIS on 25 October 2008 at 0015 UTC. Panel (a) is a false 
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color composite of one MODIS granule, showing liquid water clouds in gold, 
ice clouds in blue and white, and overlapped clouds in green. Panels (b) and (c) 
show above cloud precipitable water derived from CO 2 slicing and 0.94 pm so- 
lar absorption, respectively. Panel (d) shows the difference in the derived PW 
from these two techniques. Panels (e) shows above cloud PW derived by as- 
suming the reflecting lower cloud is located at 900 hPa and (f) the difference in 
above-cloud PW from this result and CO 2 slicing. Panels (g) and (h) show re- 
trieved values of cloud optical thickness and effective radius for liquid water 
and ice clouds. Panel (i) shows results of the precipitable water multilayer tests. 
Areas of color indicate where the tests returned a positive answer. I.e. precipi- 
table water difference was higher than 8% of the total column water vapor. 



Probability Density Function 
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a) Cloud effective radius (ice) b) Cloud optical thickness (ice) 



Fig. 2. Histograms of optical thickness and effective radius for ice clouds within the 
scene presented in Figure 1 . The effective radius histogram shows a significant 
‘shoulder’ of smaller effective radii when multilayer clouds are not removed 
from the scene. Those small effective radius retrievals come from the mixing of 
the strong liquid water cloud signal with a relatively weak cirrus signal. 
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Ice Ratio (IcR) = R 086 / R i 2 
Plant Ratio (PR) = R 0S6 / Roes 
PW fraction { PWF) = (C0 2 PW - PW 0 94 ) / TPW 
PW fraction 900 mb <PWF goo ) = (C0 2 PW - PW 0 94@900 ) / TPW 


* 

Single Layer Cloud 




i i 


Multilayer Cloud 



. 3. Flowchart for determining the presence of multilayer clouds using MODIS (col- 
lection 5). 


Height (km) 
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a) Temperature profile b) Water vapor profile 



Temperature (K) Relative Humidity (%) 

Fig. 4. Vertical profiles of (a) temperature and (b) moisture used as a basis for the for- 
ward models. The moisture profiles show saturation at levels that contain 
clouds. In this example the liquid water cloud is located at 2 km altitude and ice 
clouds are placed according to the tropopause location. 
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Fig. 5. 
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White-sky albedo as a function of wavelength for selected IGBP ecosystem 


classifications used in the forward calculations. 
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Fig. 6. MODIS multilayer cloud detection over various surfaces, water vapor content, 
and view zenith angle for a cross-section of DISORT simulations. Water vapor 
profiles are mid-latitude summer (MLS), mid-latitude winter (MLW) and mean 
tropical (TRP). Panel (a) indicates the detection of multilayer clouds over ocean 
when the lower layer liquid water cloud is placed at 2 km altitude and panel (b) 
for a lower layer liquid water cloud placed at 4 km altitude. The color bar iden- 
tifies the selection of multilayer clouds for various atmospheric profiles contain- 
ing different water vapor profiles. Panel (c) shows the MODIS multilayer cloud 
detection over various surfaces including desert, forest, and snow. Panel (d) 
shows multilayer cloud detection for various view zenith directions. 
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MODIS Collection 5 


Pavolonis-Heidinger 




profile, dark ocean, liquid cloud at 2 km 
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Fig. 7. Results of MODIS (left) and Pavolonis-Heidinger (right) multilayer cloud de- 
tection for a cross-section of DISORT simulations. Panel (a) indicates the detec- 
tion of multilayer clouds over ocean when the lower layer liquid water cloud is 
placed at 2 km altitude and panel (b) for a lower layer liquid water cloud placed 
at 4 km altitude. The color bar identifies the selection of multilayer clouds for 
various atmospheric profiles containing different water vapor profiles. Panel (c) 
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shows multilayer cloud detection for various view zenith directions. 



Fig. 8. Multilayer cloud over the western Pacific Ocean off Japan on 25 October 2008. 
Panel (a) is a true color composite of a Terra MODIS granule at 0015 UTC. 
Gray box indicates the area of interest presented in Figure 1. Panel (b) is a false 
color composite that more clearly shows the liquid water clouds (gold), ice 
clouds (blue) and multilayer clouds (green). Panel (c) is the multilayer cloud 
SDS that shows clear sky (0), single layer clouds (1), and multilayer clouds de- 
termined by various algorithm choices (2-8) (see Table 1). 
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a) MODIS Multilayer Cloud Flag (All Phases) 


b) Combined Cloud Fraction 


■ o.o 


Fig. 9. Terra MODIS monthly level-3 global products for October 2008. Panel (a) 


shows the fraction of cloudy pixels with multilayer cloud flag identified (all 
phases), and panel (b) shows the combined cloud fraction. 
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g) Cloud optical thickness 


h) Cloud effective radius 
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Fig. 10. Multilayer cloud analysis and cloud optical properties over the western Pacific 
Ocean off Japan as acquired by Terra MODIS on 25 October 2008 at 0015 
UTC. Panel (a) is a true color composite of one MODIS granule and panel (b) is 
a false color composite, showing liquid water clouds in gold, ice clouds in blue 
and white, and overlapped clouds in green. Panel (c) is another false color com- 
posite showing ice clouds in pink and liquid clouds in blue. Panels (d), (e), and 
(f) show multilayer cloud identification using the MODIS collection 5, Pavolo- 
nis-Heidinger, and Nasiri-Baum algorithms, respectively. Panels (g) and (h) 
show retrieved values of cloud optical thickness and effective radius for liquid 
water and ice clouds. 



